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Z: ■ ABSTRACT 

2 ! We investigate the Faraday rotation measure (RM) due to the intergalactic 

c/5 \ magnetic field (IGMF) through the cosmic web up to cosmological distances, 

using a model IGMF based on turbulence dynamo in the large-scale structure 

of the universe. By stacking the IGMF and gas density data up to redshift 

K*- I z = b and taking account of the redshift distribution of polarized background 

t;;;-j- \ radio sources against which the RM is measured, we simulate the sky map of the 

RM. The contribution from galaxy clusters is subtracted from the map, based on 

i„^ ^ several different criteria of X-ray brightness and temperature. Our findings are as 

O ! follows. The distribution of RM for radio sources of different redshifts shows that 

the root-mean-square (rms) value increases with redshift and saturates for -z > 1. 

The saturated value is RMj-ms ~ several rad m~^. The probability distribution 

function of |RM| follows the lognormal distribution. The power spectrum has a 



;-i ' broad plateau over the angular scale of ~ 1 — 0.1° with a peak around ~ 0.15°. 

The second-order structure function has a flat proflle in the angular separation 
of > 0.2°. Our results could provide useful insights for surveys to explore the 
IGMF with the Square Kilometer Array (SKA) and upcoming SKA pathfinders. 

Subject headings: intergalactic medium — large-scale structure of universe — 
magnetic fields — polarization 
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Introduction 



The cosmic web of filam ents and clusters of galaxies, which is predicted in the highly 

, is filled with ion ized plasma 



successful ACDM cosmo logy (iBond et a. 



galactic medium (IGM) dCen fc Ostrikerlll999[ 



1996 



i.e., the inter- 
Kanget al.ll2005f ). The hot gas with T > 10^ 



K is found mostly in the intracluster medium (ICM) and cluster outskirts, and the gas with 
10^ K < T < 10^ K, which is also referred as the Warm Hot Intergalactic Medium (WHIM), 
is distributed mostly in filaments. They were heated mostly by cosmological shock waves 
which were formed in the course of the l arge-scale structu r e (LSS) formation of the universe 



tevuet al.1 12008 : 



2008 



Vazza et al. 



Pfroinmer et al.l 120061 : iKang et al.l 120071 : ISkillman et al.l 120081 : iHoeft et al. 
I2OO9I ). The diffuse gas with T < 10^ K resides mainly in sheetlike struc- 



tures and voids. 

The IGM is expected to be permeated by magnetic fields just as the interstellar medium 
within galaxies. Theoretical studies have predicted the existence of the intergalactic magnetic 
field (IGMF). A number of mechanisms that can create seed magnetic fields have been sug- 



geste d; besides processes based on infiation and phase tran sitions (see, e.g. , IWidrow et al. 



2010l . for review) and plasma physical processes (see, e.g., IRvu et al.l I2OIOI for review) in 
the early universe, a strophysical process e s include the field g enerations duri ng the reioniza- 



tion of the unive rse (JGnedin et al.l I2OOOI : iLanger et al. 



2005), by first stars (Xu et al. 



well as the leakage of mag netic fields and cosmic ray particles from galaxies ( iDonnert et al. 



Ando et al.ll2010l ). and at cosmological shock waves (JKulsrud et al.lll997t IRvu et al.lll998l) 



2008 



as 



2OO9I : iMiniati &: Bellll201l[ ). The seed fields can be furth er amplified by fiow motions induced 



by the hierarchical clu stering during the LSS formation ( iKulsrud et al.lll997l : lRyu et al.ll2008 
Schleicher et al.lboioh . 



Faraday rotation, the rotation of the plane of linearly-polarized radio emission due 
to the birefringence of magneto-ionic media, provides an observational means to exploring 
the IGMF. There have been a number of studies of t he intracluster raagnet ic field (ICMF) 
through observations of rotation measure (RM) (see ICarilli fc Taylorl l2002l for a review). 
For instance, RMs of hundreds rad m ~^ were observed in clusters, indicating an averag e 
strength of the ICMF to be ~ 1-10 ^G Jciarke et al.lboOll : Iciarkelbooi iGovoni et aPboioh . 
In addition, RM maps of clusters were analyzed to get the power spectrum of turbulent 
magnetic fields in the ICM; for instance, a Kolmogorov-like spectrun a with a bending at a 
few kpc scale was found in the cooled core region of the Hydra cluster (jVogt fc EnBlinll2005l ). 
and spectra consistent with the Kolrnogoro v spectrum were reported in the wider ICM fo r 
the Abell 2382 cluster (JGuidetti et al.ll2008l ) and for the Coma cluster ( JBonafede et al.ll2010l ). 



Studies of RM o utside clusters, through the cosmic web, are, on the contrary, still scarce 
(e.g., IXu et al.ll2006l ). because detecting small RM is difficult with current observational fa- 
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cilities, and also removing the galactic foreground is not a trivial task. Neverthe less, recently, 



constraints for the RM through the LSS have been discuss ed by several au thors (iTavlor et al. 
2OO9I : iMao et al.1 I2OI0I : ISchnitzeleil boiol: Istil et~aD boilh. For instance, ISchnitzeleil J201oh 



argued that in the catalog of iTaylor et al.l (120091 ) . the extragalactic contribution to the width 
of RM distribution would be ctrm.eg ~ 6 rad m~^. But the measurement error is still large, 
CerrRM ~10 rad m~^; that is, so far, the nature and origin of the IGMF outside clusters 
has not been well constrained with RM studies. However, the next generation radio inter- 
ferometers including the Square Kilometer Array (SKA), and upcoming SKA pathfinders, 
the Australian SKA Pathfinder (ASKAP) and the South African Karoo Array Telescope 
(MeerKAT), as well the Low Frequency Array (LOFAR) will enable u s to investigate the 
IGMF outside clusters with high-sensitivity RM observat ions (see, e.g., ICarilli fc Rawlings 
2004J : ISecklbood : iKrause et allboooi ICaensler et aPboioh . 



Theoretical predictions for the R M due to the IGMF have been made with model IGMFs 



and s i mulations of the LSS forrnation (iRvu et al. 



2OO9I : iDolag fc Stasyszynl l2009t IStasyszyn et al. 



1998 



2010t lAkahori fc Rvul 120101 ). Especially 



Dubois fc Tevssiedl2008l : ICho fc Rvu 



based on a physically motivated model, in which a part of the gravitational energy re- 
leased during the LSS formation is transfer red to the r nagne tic field energy as a result of the 
turbulent amplification of weak seed fields, iRyu et al.l (|2008[ ) (R08) proposed that the IGMF 
follows largely the matter distribution in the cosmic web and the mean strength would be 
(B) ~ 10 nG in filaments in the local universe at redshift z = 0. Study ing various charac- 
teristic length scales of magnetic fields in MHD turbulence simulations, ICho fc Ryul (|2009[ ) 
(CR09) proposed that the RM coherence length of the IGMF in filaments would be a few 
X 100 h~^ kpc. Based on the model IGMF of R08, CR09 predicted that the-root-mean- 
square (rms) value of RM through a single filament would be of ord er ^ 1 rad m~ ^. And 



using the same model IGMF, we simulated RM in the local universe (jAkahori fc Ryull2010l ) 
(ARIO). We found that with the path length larger than the coherence length of the IGMF, 
the inducement of RM is a random walk process, but the resultant RM is dominantly con- 
tributed by the density peaks along the line of sight (LOS). The rms value of RM through a 
single filament was estimated to be ~ 1 rad m~^, which is in agreement with the prediction 
by CR09. We also found that the probability distribution function (PDF) of |RM| follows 
the log-normal distribution, and the power spectrum of the RM peaks at a scale of order 
~ 1 /i~^ Mpc. 

In this paper, we extend the RM study of ARIO in the present-day, local universe by 
including the cosmological contribution. By taking account of the redshift evolution of the 
LSS and the IGMF as well as the redshift distribution of radio sources against which RM 
is measured, we calculate the RM through the cosmic web up to z = 5. We then examine 
the spatial and redshift distributions of the RM and their statistics. In Sections 2 and 3, 
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we describe our models results, respectively. Discussion is in Section 4, and summary and 
conclusion follow in Section 5. 



2. Models 

2.1. Large-Scale Structure of the Universe 

For the LSS of the universe, we used structure formation simulations of a ACDM universe 
with the following values of cosmological parameters: flho = 0.043, fimo = 0.27, Qao = 0.73, 
h = Ho/{100 km/s/Mpc) = 0.7, n = 1, and Ug = 0.8. The sir nulations were p erformed 



using a particle- mesh/Eulerian, cosmological hydrodynamic code ( iRyu et al.lll993l ). A cubic 
region of comoving volume (100 h^^Mpc)^ was reproduced with 512^ uniform grid zones for 
gas and gravity and 256'^ particles for dark matter, so the spatial resolution is 195 h"^ kpc. 
Sixteen simulations with different realizations of initial condition were used to compensate 
cosmic variance. They are the same set of simulations used in ARIO. 



2.2. Intergalactic Magnetic Field 

As in ARIO, we employed the model described in R08. It assumes that turbulent-flow 
motions are induced via the cascade of the vorticity generated due to cosmological shocks 
during the formation of LSS, and the IGMF is produced as a consequence of the amplification 
of weak seed fields of any origin through the stretching of field lines by the flow motions. 
Then, it can be modeled that a fraction of the turbulent-fiow energy, eturb, is converted to 
the magnetic energy, Eb, as 

^B = ( I e^turb, (1) 

V teddy/ 

where the conversion factor, 0, depends only on the eddy turnover number, t/teddy Here, 
the eddy turnover time is defined as the reciprocal of the vorticity at driving scales, teddy = 
1/i^driving (a; = V X -u). The local vorticity and turbulent-flow energy density were calcu- 
lated from the data of the structure formation simulations described above, and the age of 
the universe at the redshift z was used for t. The functional form for the conversion factor 
was derived from a separate, incompressible, magnetohydrodynamic (MHD) simulation of 
turbulence dynamo (ROB). Then, the magnetic energy density was calculated according to 
Equation [H and the strength of the IGMF as B = v^SttFb- The average strength of the 
resuhing model IGMF for the WHIM is {B) ~ 10 nG or {pB)/{p) ~ 0.1 /xG at z = 0. For 
the direction of the IGMF, we used that of the passive fields from structure formation simu- 
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l ations, in which weak seed fields were generated through the Biermann battery mechanism 



( iBiermannl Il950[ ) at cosr aological shocks and evolved passivel y, ignoring the back-reaction 



along with flow motions (JKulsrud et al.l 119971 : iRyu et al.lll998l ) 



The validity of our model IGMF was discussed in details in ARIO. It was argued that 
our model IGMF would produce reasonable results for the RM through the cosmic web in the 
local universe. The simulated passive flelds reproduce the directions that show the expected 
correlation with those of vorticity. The RM coherence length for the IGMF in fllaments is 
estimated to be a few to several xlOO kpc at z = 0, which agrees with the estimation of 
CR09. It is a few times larger than the grid resolution of our simulations, 195 h~^ kpc. 
On the other hand, the coherence length for the IGMF in clusters is expected to be a few 
xlO kpc (CR09), which is smaller than the grid resolution. So the application of our model 
IGMF to clusters would results in an erroneous estimation of RM. However, here we study 
the RM through the cosmic web outside clusters, and so we will remove the contribution 
from clusters to RM (see Section 2.8). 

Our model predicts that the IGMF for the WHIM was a bit stronger in the past; for 
instance, (B) ~ 30 nG for the gas with 10^ K < T < 10^ K at ^ = 5 (see Figure 2 of R08). 
It is because the density of the WHIM was higher in the past, although the magnitudes 
of vorticity and the vortical component of velocity were smaller. Our model also predicts 
stronger IGMF for the hot gas with T > 10^ K in the past. However, the IGMF averaged 
over the entire computational volume was weaker in the past, because the volume and mass 
fractions of strong-fleld regions were smaller. 



2.3. Data Stacking 

To reproduce the cosmic space, we stacked simulation boxes u p to 2: = 5, basicall y 



following the conventional manner of cosmological data stacking (e.g., Ida Silva et al.ll2000l ). 
We took outputs at Zout = 0, 0.2, 0.5, 1.0, 1.5, 2.0, 3.0, and 5.0, and put the outputs of 
the closest redshift for stacking boxes; 56 boxes were required to reach z = 5. The stacking 
boxes were randomly selected from sixteen simulations and randomly rotated to avoid any 
artiflcial coherent structure along the LOS. If a LOS went out of the boundary of a box, 
we applied the periodic boundary condition in which we replicated the box across the LOS. 
We also carried out the calculation with the open boundary condition in which we did not 
replicate the box but put different one from sixteen simulation outputs, and conflrmed that 
the statistical properties of RM we examined are not sensitive to the boundary conditions. 
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2.4. Observer Locations 



Our Galaxy is located in the Local Group. Since the Local Group is probably filled with 
the magnetized WHIM, the observed RM could contain a contribution from the IGMF of 
the Local Group (see Section 3.1). It would have been ideal to place the "observer" where 
the IGMF and gas density are similar to those in the Local Group. Unfortunatel y, however 



not m uch is known about the physical states of the Local Group. Hence, following iDas et al. 



( I2OO8I ). instead, we selected groups of galaxies identified in the simulation data that have 
similar halo gas temperatures to the Local Group. We considered two cases: groups at 
z = with 0.05 keV < kTx < 0.1 5 keV and groups at z = with 0.05 keV < kTx < 0.5 
keV (JRasmussen fc Pedersenll2001l ). Here, Tx is the X-ray temperature (see Section 2.7). 
Assuming that these groups are located in a physical environment similar to that of the 
Local Group, we placed "mock observers" at the center of the groups. 



2.5. Survey Setup 

We aim to simulate wide-field, high- sensitivity RM surveys that will be achieved with 
future radio interferometers. For instance, the upcoming POSSUM (Polarization Sky Survey 
of the Universe's Magnetism) project, one of the major surveys planned for the ASKAP, will 
produce an RM map that will be a dramatic improvement over existing maps. The RM map 
from the POSSUM will have a field-of-view (FOV) of 6*^ ~ 30 deg^ and an average source 
separatio n of (A^) :^ 0.1° (~ 6 arcmin) corresponding to an RM grid of ~ 100 RM deg~^ 
(see, e.g.. lGaensler et al.ll2010l ). The SKA w ill survey ^ 10^ polarized radio sources over the 



whole sky with (A^) ~ 1 arcmin (see, e.g. 

map of RM, we adopted 9^ = 200 deg^ = (14. 14)^ deg^ and A9 = ^^200 degV20482 



Beclg |2009|) . Intending to produce a simulated 

0.414 



arcmin with one source in each of 2048^ equally spaced pixels. Specific simulations for the 
SKA survey and the ASKAP POSSUM that consider the feasibility of observations such as 
the effects of the measurement error, coarser spacing, randomly placed sources, intrinsic RM, 
and so on, as well as the galactic foreground will be discussed in separate papers. 

For the redshift distribution of radio sources, we employed a d istribution based on the 
estimation of the detectable radio galaxies by the SKA and ASKAP (jWillman et al.ll2008L see 
also Figured]), in which the observed FR I and FR II galaxies are taken into consideration. 
The number of available sources drops substantially at 2; > 5, and that is the reason why 
the data stacking was made up to z = 5 (see Section 2.3). 
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2.6. Faraday Rotation Measure 



RM is a measure of tlie Faraday rotation of polarized radio emission against a back- 
ground source, defined by 



^ e^ r= „ A(/) 



™^^=2^^^^^ -ei^ii^rf/, (2) 

^obs ^vrmgC Jo ^obs 

where ip is the rotated angle of the plane of polarized radio emission, Ig is the path length up 

to the source, Ug is the thermal electron density, B\\ is the LOS component of the magnetic 

field, Aobs is the observed wavel ength, and A(/) is the way elength along the path length 



of polarized radio emission (e.g., iRybicki fc Lightmanlll979[ ). For a source at cosmological 
distance of redshift Zs, it is modified to 

RM= -^ /' ' {l + z)-'n,{z)B\\{z)dl{z), (3) 

where ne{z), B\\{z), and l{z) are the quantities in proper coordinates. The path length is 
given as 

dl{z) = / "^^ (4) 



for the flat universe with fi^o + ^ao = 1 (e-g-, IPeebled 119931 ). 



We calculated Equation ([3]) by counting the contribution from computational grid zones 
along stacked boxes up to mock sources as 

RM (rad m'^) = 8.12 x 10^ ^ {^ + ^i)~^ " ^e(^i) ■ B\\{zi) ■ Al{zi), (5) 

where Ns{zs) is the number of grid zones up to sources at Zg- The density and magnetic field 
strength are in units of cm~^ and /iG, respectively. For Al{z), we used the proper size of 
grid zone, which is 0.195 h^^{l + z)''^ Mpc. With the simulation FOV of 6"^ = (14.14)^ deg^, 
the maximum tilt of LOSs relative to the coordinates of simulation box is cos(6'/2) = 0.99, 
so we ignored the effect on Al{z). 



2.7. X-ray Emission 

In the hot gas with T > 10^ K, X-ray emission is produced mainly by thermal bremsstrah- 
lung. We hence considered only the bremsstrahlung emission from electrons, and neglected 
line emissions from ions. The emissivity of thermal bremsstrahlung can be expressed as 

£«f (erg S-' cm-3) = 7.77 x IQ-'^^T^'^^nl ■ H g{T, u) exp (-^) dp, (6) 



where T, Ug, and u are in units of K, cm ^, and Hz, respectively, and h p and k^ are the 



Planck and Boltzmann constants, respectively ( iRybicki &: Lightniadll979l ). We adopted the 
approximate Gaunt factor, g ^ .dlhpU / k^T)~^'^ , and used the bolometric emissivity, that 
is, 1^1 = and z/2 = 00. 

For the X-ray temperature of clusters and groups, we calculated the X-ray emissivity- 
weighted temperature as 

Tx= I e^TdV f e^dV, (7) 

over a spherical volume of comoving radius 0.5h^^ Mpc. 

For the X-ray surface brightness and surface temperature up to redshift z, we first sought 
all the grid zones, k, and their proper volume, \4, which enters within the angular beam of 
(A^)^. Note that A6 = 0.414 arcmin corresponds to the comoving size of a grid zone 195 h~^ 
kpc at 2; ~ 0.6; hence perpendicular to a LOS, less then one zone enters at lower redshifts, 
but more than one at higher redshifts within (A6')^. For each grid at redshift z^, the X-ray 
luminosity was calculated as 

Lk= f e^'dV. (8) 

Then, the X-ray surface brightness was calculated as 

up to 2 

^ (Aey ^ ^T,diz^)\\ + z^r ^ ^ 

and the X-ray emissivity-weighted surface temperature was calculated as 

up to ^ rji J 

T* = V ^^ '' (10) 

^ S*^{Aey ^ A-nd^izkYil + ZkY' ^ ^ 

where Txk is the temperature of grid zone k. Here, dc{zk) is the comoving distance up to z^, 

which is given as 

,/ X c /■"''= dz , , 

^0 Jo Vl^mo(l + ZY + ^Ao 



2.8. Subtraction of Cluster Contribution 

The contribution from clusters was subtracted from calculated RM based on different 
criteria using X-ray brightness and temperature. The simplest approach would be to exclude 
the volume associated with clusters in the cosmic space. We employed the following two 
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criteria. In the first model, labeled "TM7" , all the hot gas with T > 10^ K was excluded. In 
the second model, labeled "CLS", the spherical region of comoving radius lh~^ Mpc around 
the clusters and groups with Tx > 2 keV was excluded. 

For simulations of RM surveys, however, a better approach would be the exclusion of 
the pixels that include clusters in the projected sky map. We would like to subtract the 
contribution from all clusters, but it is still hard to detect faint X-ray clusters. So we tried 
the following criteria. In the model labeled "TS8", all the pixels with TJ > 10^ K and 
Sx > 10~^ erg s~^ cm~^ sr~^ were exc luded; 10~^ erg s~^ c m~^ sr~^ is close to the detection 



limit of current X-ray facilities (e.g., iHoshino et al.ll2010[ ). And in another model labeled 
"TSO", the pixels with T|- > 10^ K and S^ > 10"^° erg s"^ cm'^ sr'^ were excluded; 10"^° 
erg s~^ cm~^ sr~^ is intended to mimic the improved detection limit of future X-ray facilities. 
Here, Sx and TJ up to z = 5 were used. Finally, we also considered the case where no volume 
or pixel was excluded, labeled "ALL", for comparison. 



3. Results 

3.1. Contribution from Groups around Observers 

We first examine the contribution to RM from the local IGMF around observers. For 
the two temperature range of groups containing observers, 0.05 keV < kTx < 0.15 and 0.05 
keV < kTx < 0.5 keV (see Section 2.4), we set up 200 mock observers and calculated RM 
by integrating along LOSs up to ~ 0.5 h~^ Mpc (21/2 grid zones) from observers. Figure 
[2] plots the PDF of |RM|, the absolute value of resulting RM. The PDF peaks around 
~ 2 X 10^3 rad m'^ and ~ 6 x lO'^ rad m'^ for 0.05 keV < kTx < 0.15 keV and 0.05 
keV < kTx < 0.5 keV, respectively. The rms values are ~ 5 x 10^^ rad m^^ and ~ 6 x 
10^^ rad m~^, respectively. The peak values can be estimated roug hly as follows. Our model 



IGMF predicts Epeak ~ 10"^ /xG for groups fsee iDas et al.ll2008h . With n^ ~ 10"^ - 10"^ 
cm~^ for groups, and the coherence length of magnetic fields comparable to the size of grid 
zones, i.e., /coherence ~ 195 h^^ kpc, we get RMpeak ~ 10^'^ — 10~^ rad m~^. This may 
indicate that the contribution from the local IGMF is affected by the grid resolution. But 
as we will see in Section 3.3, this contribution makes only a small, negligible fraction of the 
RM through the cosmic web. Below we show only the results with 0.05 keV < kTx < 0.15; 
those with 0.05 keV < kTx < 0.5 keV are statistically indistinguishable. 
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3.2. Two-Dimensional Map 

We produced the simulated sky map of RM (using Equation ([5])) with different cluster 
subtraction models, as well as the accompanying map of 5*^ (using Equation Q) and TJ 
(using Equation (ITO|) ). for 200 different mock observers with different stacks, each of which 
has 2048^ pixels in a FOV of 200 deg^, as described in Section 2.5. The statistics below were 
obtained with the maps. In this section, we show maps for a typical case. 

Figure [3] shows RM maps for ALL and cluster subtraction models CLS and TS8, and 
Sx and TJ maps. To see the dependence on the redshift depth, the maps integrated up 
to five different epochs are presented. With TS8, about ~ 15% of pixels were subtracted. 
At the lowest redshift depth. X-ray bright sources, which can be identified as clusters or 
groups, are apparent. They have the angular size of up to a couple of degrees and the peak 
value of |RM| > 100 rad m~^. With CLS, some, but not all, contribution to RM from those 
clusters and groups was excluded. With TS8, most contribution from clusters and groups 
was removed. On the other hand, the contribution from filamentary structures was kept; 
filamentary structures with |RM| ~ 0.1 — a few rad m~^ are seen in CLS and TS8 as well as 
in ALL. At higher redshift depths, the value of |RM| is larger. This is attributed to the fact 
that the inducement of RM is a random walk process (ARIO); the value increases with the 
increasing path length. Also at higher redshift depths, the angular size of RM structures is 
smaller, since the distance is larger and the proper size of filaments is smaller. At the depth 
of 2 = 5, most of nearby filamentary structures are smeared out, and the sky distribution of 
RM appears close to be random. 



3.3. PDF and Root Mean Square 

We quantified the RM in 2048^ pixels of 200 maps, so in total 2048^ x 200 pixels, first 
with the PDF and rms value. Figure H] shows the PDF of |RM| for three different redshift 
depths, z = 0.05, 0.3, and 5, for ALL, CLS, and TS8. As for the |RM| of the local universe 
(ARIO), the PDF follows the lognormal distribution. For TS8, the peak value shifts from a 
few xlO~^ rad m~^ at the depth oi z = 0.05 to ~ 1 rad m~^ at z = 5. Compared to CLS 
and TS8, ALL has a high-value tail due to the contribution from clusters and groups. 

Figure [5] shows RMrms, the rms value of RM, as a function of the redshift depth for 
ALL and all cluster subtraction models. RM^ms increases steeply for 2; < 1, and saturates 
for z > 1. The saturated value of RMj-ms is ~ 10 rad m~^ for CLS and ~ 7 — 8 rad m~^ 
for TM7, TS8, and TSO, while it is ~ 40 rad m~^ for ALL. Again the high value for ALL is 
attributed to clusters and groups. So if we observe the RM through the cosmic web outside 
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of clusters, we would get ~ 7 — 8 rad m ^ for RMrms- 

We note that the saturated value depends on the redshift distribution of radio sources, 
but only weakly. In our model the redshift distribution peaks at 2; < 1 (see Figured]). We 
carried out comparison runs, in which we put all radio sources at z = 5. Then, the saturated 
value of RMrms increased slightly; specifically, for TM7, the value increased from ~ 8 to ~ 9 
rad m~^. 

The saturation for 2; > 1 is mostly due to the convergence of the path length at 2; ~ 1 
and the (1 + z)~^ dependence of KM (see Equation ([5])). To see this, we calculated the path 
length and also the column density across the WHIM with 10^ K < T < 10^ K as a function 
of the redshift depth. Figure |6] shows the resulting path length and column density. The 
path length saturates at ~ 100 Mpc for z > 1. On the other hand, up to z = 1, the column 
density reaches ~ 40 % of the saturated value. But the contribution to RM from higher 
redshift should diminish due to the (1 + z)~'^ dependence. 

We can roughly estimate the saturated value of RMrms, as follows. Since the inducement 
of RM is a random walk process, we may take the value as ~ RMrms, filament x (^filaments) ^, 
where RMrms, filament is the rms RM induced by a single filament and A^aiaments is the number of 
encounters of filaments. ARIO estimated that R Mrms, filament ^ 1 - 5 rad m~^. Since the typical 



width of filaments is several h~^ Mpc (see, e.g.. lRyu et al.ll2003l : lKang et al.ll2005l ). the path 
length of ~ 100 Mpc corresponds to Armament ~ 25. Then, we get RMrms ~ 1-5 x a/25 ~ 7.5 
rad m~^, which agrees well with those for TM7, TS8, and TSO. 



3.4. Power Spectrum and Structure Function 

We also calculated the two-dimensional power spectrum with the sky map of RM. Figure 
[7] shows the resulting spectrum, averaged over those from 200 maps, for redshift depths 
z = 0.05, 0.3, and 5. For our best models TS8 and TSO, at the depth oi z = 0.05, the power 
spectrum peaks at ~ 1°, which corresponds to the proper length of ~ 2.5 h^^ Mpc at that 
redshift. As the redshift depth increases, the power adds up mostly at smaller scales. At 
z = 5, the power spectrum shows a broad plateau over ~ 1 —0.1° with peak around ~ 0.15°. 
The peak angle corresponds to the proper length of ~ 3 h^^ Mpc at z = 1 and ~ 2.5 h^^ 
Mpc at z = 5. With the typical radius of filaments of a few h~^ Mpc, we argue that the 
peak in the power spectrum reflects the scale of filaments. 

The point that the power spectrum peaks at a small angular scale of ~ 0.15° could 
be the key to the search of the RM due to the IGMF in future RM surveys such as the 
SKA survey and the ASKAP POSSUM; it can be used to remove the galactic foreground. 
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Towards the galactic poles, the galactic foreground is e xpected to have t he peak of the power 



spectrum at a much larger scale of ~ 10° or so (e.g., iFrick et al.ll200ll ). Then, it would be 
possible to separate the extragalactic component from the galactic foreground in the Fourier 
space, when RM surveys become available. The issue of the galactic foreground, along with 
other uncertainties, in the study of the extragalactic component of RM will be discussed in 
separate papers. 

From the observer's point of view, the structure function (SF) is a statistical quantity 
which is easier to be obtained. The the n-th order SF is defined as 



Snir) = {\RMix + r)-RM{x)\'')s 



(12) 



with r = |r|, where the subscript indicates the averaging over the data domain. It can be 
easily calculated for the data of arbitrary domain with irregular sampling intervals. Figure [S] 
shows the second-order SF (n = 2) for the depth of z = 5 in our model, again averaged over 
those from 200 maps. The SF increases slightly at small angular separations, < 0.2°, and 
saturates and stays flat at larger separations. The second-order SF is basically the Fourier 
transform of the power spectrum. The saturation at ~ 0.2° reflects the fact that the power 
spectrum has the peak around that angular scale, which is indeed the case as shown in Figure 

m 



Recently, for instance. iMao et al.l (J2010|) and lStil et al.l (1201 if ) reported the second-order 
SF from RM surveys towards the galactic poles; their data are also shown in Figure [HI The 
SF towards the south pole is relatively fiat on the angular scale of ~ 1 — 10°, while it decreases 
with decreasing angular separation at smaller scales. On the other hand, the SF towards 
the north pole decreases with decreasing angular separation on the scale of ~ 1 — 10°, while 
it is very noisy but looks fiat at smaller scales (the noise is partly attributed to the Coma 
cluster in the field ). The SF toward s the south pol e is so mewhat larger than that toward 
the north pole. As ISchnitzelen ( l2010l ) and IStil et al.l ( 1201 ll ) noted, however, the error in the 
observed SFs is still too large to allow any concrete conclusion. Nevertheless, our results 
suggest that a substantial fraction of the RM towards the galactic poles may be attributed 
to the IGMF. As a matter of fact, our estimation of the extragalactic contribut ion to RM, 
~ 7 — 8 rad m~^, agrees with the estimation of ~ 6 rad m~^ from a survey by ISchnitzeler 



mm 



4. Summary 

Measuring Faraday RM is one of a few available methods to investigate the IGMF in 
the LSS of the universe, especially in filaments, which still remains largely unknown. Using 
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a model IGMF derived from cosmological structure formation simulations, we examined the 
RM through the cosmic web outside of clusters, by stacking simulation data up to redshift 
z = 5 and taking account of the redshift distribution of radio sources. We found that the 
magnitude of the RM increases with the path length; the path length and so the value of 
RMrms increases with the redshift depth of radio sources for 2; < 1 and saturates for 2; > 1. 
The saturated value is predicted to be RM^ms ~ 7 — 8 rad m~^. We also found that the PDF 
of |RM| follows the lognormal distribution. The power spectrum of the RM shows a broad 
plateau over the angular scale of ~ 1 — 0.1° and peaks at ~ 0.15°. The second-order SF 
shows a flat profile in the angular separation of > 0.2°. 

Our results could provide useful insights for planning RM surveys with the SKA, and 
upcoming SKA pathfinders, ASKAP, MEERKAT, and LOFAR. Also we point that our 
results could be tested with such surveys. For instance, our power spectrum of RM peaks at 
the scale of ~ 0.15°. A very fine RM map with the average separation of sources ~ 1 arcmin, 
which is expected to be achieved with the SKA, should have enough angular resolution to 
see the power spectrum peak we predict. 

We thank J. M. Stil for proving the data used in Figure [8] and H. Kang for comments 
on the manuscript. The work was supported by the National Research Foundation of Korea 
through grant 2007-0093860. 
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Fig. 1. — Redshift probability distribution function of rad io sources based on the expected 
number of observable FR I + FR II galaxies by the SKA (jWillman et al.ll2008l ). 
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Fig. 2. — PDF of |RM|, the absolute value of RM, due to the local IGMF in groups where 
observers are located. Solid and dotted lines show the results with the temperature criteria 
for the groups, 0.05 keV < kTx < 0.15 keV and 0.05 keV < kTx < 0.5 keV, respectively. 
200 simulations are carried out for each case. 
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Fig. 3. — Maps of RM for ALL, CLS, and TS8 models, X-ray surface brightness, Sx, and 
X-ray emissivity- weighted average temperature, T^, integrated from ^ = up to the redshift 
depth shown. The angular size of each map is 14.14 x 14.14 deg^. 
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Fig. 4. — PDF of |RM| integrated up to the redshift depth shown. The average for 2048^ x 200 
pixels is shown. Sohd, dashed, and dotted hnes show the results for ALL, CLS, and TS8 
models, respectively. 
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Fig. 5. — RMrms, the rms value of RM, integrated up to the redshift depth z = 5, for ALL, 
CLS, TM7, TS8, and TSO models. The average for 2048^ x 200 pixels is shown. 
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Fig. 6. — Path length (dashed) and column density (solid) across the WHIM with the 
temperature of 10^ K < T < 10^ K, as a function of the redshift depth. The average for 
2048^ X 200 pixels is shown. 
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Fig. 7. — Two-dimensional power spectrum of RM for ALL, CLS, TM7, TS8, and TSO 
models. The average for 200 stacking simulations is shown. For each model, the power 
spectra integrated up to the redshift depth z= 0.05, 0.3, and 5.0 are shown. 
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Fig. 8. — Two-dimensional, second-order structure function of RM for CLS, TM7, TS8, 
and TSO models. The average for 20 stacking simula tions is show n. Also shown a re the 
second-order structure functions from Mao et al.l (I2OIOI ) (circles) and lStil et al.l ( 1201 ll ) (solid 
lines). Filled circles and thick lines are toward the south galactic pole, and open circles and 
thin lines are to ward the north galactic pole. They are the same as those in Figure 14 of 
Stil et all f l201ll ). 



